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We present a class of non-Gaussian two-mode continuous variable states for which the separability 
criterion for Gaussian states can be employed to detect whether they are separable or not. These 
states reduce to the two-mode Gaussian states as a special case. 
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PACS numbers: 03.67.-a 

Whether a quantum state is entangled or not repre- 
sents a very important question in quantum-information 
theory. Such knowledge reveals whether one can take ad- 
vantage of the non-local properties of the state in quan- 
tum protocols such as quantum teleportation [l[ and 
quantum cryptography [2j. This issue has been dealt 
with by many authors in recent years primarily in qubit 
systems where the Peres- Horodccki partial transpose sep- 
arability condition was the first method to figure out if 
a two-qubit state was separable Q. In general, for N 
qubits the solution is not known. Continuous variable 
systems have proven to be an extremely useful setting for 
quantum cryptography and communication [ij]. In these 
protocols entangled states are required and the question 
of separability arises naturally. For two-mode systems 
separability criteria for Gaussian states were established 
in Refs. [5|, [6( which proved to be both necessary and 
sufficient. More recently necessary and sufficient condi- 
tions for the partial transposition of bipartite harmonic 
quantum states to be nonnegative are formulated in Ref. 
0, [§[ and separability criteria based on uncertainty rela- 
tions for two-mode representations of SU(2) and SU(1,1) 
algebras in Refs. 0, |l(| EH can be obtained from from 
the former result as special cases. These criteria have 
particularly targeted uncovering whether non-Gaussian 
states are separable or not as previous criteria fail to de- 
tect relatively simple entangled states. 

In this paper we approach the problem from a different 
perspective. The Wigner function of the reduced state of 
a two-mode quantum state is shown in Fig. [TJ Clearly 
the Wigner function is non-Gaussian and hence the two- 
mode state from which it is derived is non-Gaussian as 
well. Say an experimentalist measures (that is recon- 
structs from a tomographically complete measurement) 
this state in the laboratory. What exactly can she say 
about the separability of the state? At first sight not 
a great deal. However we will show that in fact this 
state's separability is completely known and understood 
in terms of its interpretation as a Gaussian state. Fur- 
thermore we introduce a whole class of these states for 
which the usual two-mode Gaussian states are a special 
case. 

Gaussian states of a two-mode continuous variable 
quantum system have been much studied in the liter- 



ature and their entanglement properties are quite well 
established [H], Il3j . Here we briefly review their char- 
acterization. Let a\a be the bosonic creation and an- 
nihilation operators acting on a Fock Hilbert space as 
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the Weyl-Heisenberg commutation relations 
We can introduce position and momentum operators x = 
x* + x, p — i(a* — a) (setting h — I) which define phase 
space through the continuous range of their eigenvalues. 
For two modes these position and momentum operators 
are defined in each Hilbert space. We gather these oper- 
ators into a single vector x = (xi,pi,X2,p2) for clarity. 
By definition the Wigner function of a Gaussian state, p, 
takes the form W{x) = exp[— x'ax' T /2}/n\ / det a, where 
x' = x — (x) and ajk = ({xjXk + Xk£j))/2 - (xj){xk) 
is called the covariance matrix with (O) = Tr[pO]. The 
Wigner function thus depends only on the first and sec- 
ond moments of the position and momentum operators 
[bH ]. Furthermore, via local operations the covariance 
matrix can be brought to the simpler form 
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where bi,Ci and di satisfy 



and I ci 



c 2 |= v/(&i 1) ~ v(^2 - 1)(^2 - 1) as shown in 

Ref. [a]. The necessary and sufficient condition for the 
state to be separable is then 
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with u = qoXi 
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p2, and ql = 

l)/(bi — 1). Actually this expression simplifies to 

when both b\ — 
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Cj | otherwise. In the following section we im- 



pose an additional criterion: we will relate the degree of 
entanglement of Gaussian states to an energy of these 
states. 

Gaussian states are an important class of states, but 
some quantum protocols require [l5| the usage of non- 
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Gaussian states. These states are more difficult to handle 
mathematically than Gaussian states which only require 
knowledge of the first and second moments of a finite 
set of observables. One of the reasons why it is more 
difficult to analyse non-Gaussian states is that they are 
characterized by an infinite set of non-zero cumulants 
i.e. higher-order moments of system observables cannot 
be expressed in terms of the first and second order mo- 
ments. What is important to note is that in order to de- 
fine Gaussian states one needs to be able to construct ob- 
servables from creation and annihilation operators which 
satisfy the Weyl-Heisenberg commutation relations. It 
is therefore possible to choose other operators satisfying 
these commutation relations and use these to construct 
new observables whose eigenvalues define a completely 
different phase space. In particular we can choose multi- 
photon operators like 
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as in Ref. pi], satisfying [iW , iWt] = II where 
N = fita, [[N/k]] = EJ[n/fc]]|n)<n|, (N - k)\/N\ = 
Y^ n fo~ fc ^ ! | n){n\, and [[n/Jfc]] denotes the largest posi- 
tive integer less than or equal to n/k . They are con- 
structed in such a way as to create or annihilate k 
photons at a time. These operators can also be inter- 
preted as acting on a "multi-photon Fock space" , TL, as 
i«t|„) fe = ^+T\n+l) k , A( k 1\n) k = Vn|n-l) fc 
where the subscript on the multi-photon number state 
| n) fe indicates we are referring to those states satisfying 
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Their action on the usual Fock space is given by 

nk + m) = yn + 1| (n + l)k + m) 
A (k) \nk + m) = y/ri\ (n - l)k + m), 



(5) 



with n > and < m < k—1. Thus, for the multiphoton 
number operators — A^ A^ each eigenvalue is k- 
times degenerated (including the vacua, i.e. i( fe ) | m) = 
for m < k). 

As before we can construct position and momen- 
tum operators (acting on H. only) for the two modes 
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A^ ) and the eigenvalues of these op- 



i(fe) 



erators define the phase space. In this phase space we 
have Gaussian states whose Wigner functions are Gaus- 
sian while if we represent the states using a Wigner func- 
tion in the usual phase space the states are highly non- 
Gaussian. The separability criterion from Ref. [6[ carries 
over here so that given a state has a covariance matrix 
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we know it is separable if 
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otherwise. Of course the local operations required to 
transform the state under consideration to one which 
can be represented by such a covariance matrix are non- 
Gaussian operations in terms of the creation and annihi- 
lation operators a\, a\, eij, a%. Experimentally the opera- 
tions could be constructed as proposed in Ref. [ItJ ■ 

As an example we can take the multi-photon two- 
mode squeezed vacuum state which has the form | VP) = 
yl — 7 2 J2 n 7™ I n ) fe l n ) fe - The re duced state in either 
mode is given by p^ = (1 — 7 2 ) ^2 n J 2n \ n)k{n |, a ther- 
mal state. The Wigner function of this state in the 
(X^ , p[ k ^ ) reduced phase space is shown in Fig. [2] and 
has the expected Gaussian form. However a simple cal- 
culation shows us that 
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giving a reduced density matrix p = (1 — 
7 2 ) E n 7 2 "| fcn) (fcn |. The Wigner function for this 
state in the (xi,pi) reduced phase space is plotted in 
Fig. [1] showing that it is clearly non-Gaussian and even 
contains negative parts. 

The covariance matrix of the two-mode multi-photon 

squeezed state is as in Eq. ([6]) with B^ = = 
cosh2r and C± = sinh2r = —C% where 7 = tanhr. 
The separability criterion in Eq. ([7]) reads exp(— 2r) > 1 
and the state is entangled for r > 0. How then does 
the separability criterion relate to measurements of ob- 
servables in the usual Fock space? For example we need 



to know what B 



(fe) 



is and to this end we must measure 
(x[ k ^ 2 ) and (x[ k ^). For the case k = 2 we find these 
expectation values in terms of the operators a\ , a\ are 
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and 

(ii) 

with N± = a\d\. While we readily concede that it is a 
non-trivial task to experimentally measure these expec- 
tation values we are motivated by what one is able to say 
about the separability of a state when full tomography 
has been carried out and the result is a state such as that 
in Fig. 1. 

For a given squeezing parameter 7 the two-mode 
squeezed vacuum state | $) = yl — j 2 J^ Tl 7"] n, n) and 

its multi-photon equivalent | = \J\ — 7 2 J2 n 7™ I n > n ) k 
will posess the same degree of entanglement, a fact 
most easily seen using the the von Neumann entropy of 
the reduced state of either mode, S = — Trpilnpi = 
— Trpi l n Pi = S( k \ It should be noted however that 
uncovering the two-mode entanglement present in each 
of the above states requires measurement of different ob- 
servables. In what follows we will compare two states of a 
two- mode system at a fixed energy value [l3T ] . In terms of 
the average energy (E) = ($ | J^i \ ^) w ith = 
the multi-photon squeezed vacuum state with the same 
degree of entanglement has average energy (E^) = 
k(E). If we fix the energy of both states at (E) = 2h 
where n = ($ |iVj| <fr) then the reduced state of | $) can 
be written as p\ = J2 n Pn\ n ) ( n I withp„ = n n /(l+n) 1+n 
and the reduced state of | ^) can be similarly written as 

P [ i ] = Y,nPn ] \n) k {n\ withpW = (|)"/(1 + 
Thus it is clear that for fixed energy the usual two-mode 
squeezed vacuum is more entangled than its multi-photon 
counterpart. Intuitively this makes sense given that the 
nature of the multi-photon state does not allow certain 
quantum correlations to exist; for instance, upon mea- 
surement of the state in the joint number basis there is 
zero probability to obtain the result | kn + m, kn + m) 
for m < k and n E Z + . This result also ties in with 
the fact that among all continuous variable states with 
a given fixed energy, the maximally entangled states are 
Gaussian ( 181 ] - 

For mixed states we compare a fc-photon two-mode 
mixed state with a usual two-mode mixed state having 
the same average energy (E). To get an intuitive sense 
of the difference between two such states we look for the 
minimum purity (defined as P(u>) = Tilu 2 ) allowed for 
the fc-photon mixed state given this energy (E). The 
dependence on k is 



corresponding to a tensor product of two multi-photon 
thermal states. As fc increases the minimum purity in- 
creases asymptotically toward 1 so that the fc-photon 
mixed states tend toward the vacuum in the limit fc — > 00. 
To re-enforce this point in Fig. [3] we plot the maximally 



entangled multi-photon Gaussian mixed states, see Ref. 
[l3l ] . for various values of fc, all at a fixed mean energy 
(E) = 1. Thus we can say that the maximally entangled 
fc-photon mixed states are less entangled than those for 
fc = 1. For general mixed states this statement is not 
always true. 

We have presented a large class of non-Gaussian states 
for which the existing separability criterion for Gaussian 
states can be employed in order to detect their entan- 
glement. In order to clarify our results we recall that 
an arbitrary unitary transformation U : H — > H result- 
ing in an "annihilation" operator bjj = UaU' can be 
exploited to define another class of non-Gaussian states 
Cu = {q ■ Q = U*qU,q € Q}, where Q denotes the set 
of "standard" Gaussian states. Due to the fact that uni- 
tary transformations preserve an operator's spectra and 
the commutation relations, the operators bu,b\j form a 
representation of the Weyl-Heisenberg group and Gaus- 
sian states with respect to these operators can be defined. 
In fact, as it was pointed out in Ref. Q the inequality in 
Eq.([2]) provides a sufficient separability criterion for all 
operators jx 1 v that are locally unitary equivalent to u, v, 
i.e. jj, = UuU', v = UvU\ respectively. Moreover, these 
inequalities provide the necessary conditions for entan- 
glement for all Gaussian states defined with respect to 
new phase-space coordinates, i.e. for all qeC\j. 

We have to note that even though the multi-photon 
non-Gaussian states analyzed in our paper seem to be 
of the similar form as discussed above, there is a signif- 
icant difference. The operators a and A k are not mu- 
tually related by a unitary transformation in the above 
sense (for more details see Ref. .19]). In fact, the num- 
ber operators Nk and n have different spectra (Nk is 
degenerated). The construction in our case is based on 
the fact that the (semi)infinitc Hilbert space of the orig- 
inal harmonic oscillator can be expressed as a finite di- 
rect sum of (semi)infinitc Hilbert spaces that are isomor- 
phic to the original one, i.e. 7i — Hq © • • • © Hk-i- 
Here Hj is a linear span of vectors H 3 \nk + j) = 
\n)j,k S Tij (n — 0, ...,00). Physically this means 
that we are restricted to states belonging to the sub- 
space spanned on photon number states separated by 
a fixed energy ktku (representing the energy of fc pho- 
tons). The vacuum for Tij is represented by the state 
I + j) = \j) G H The linear spaces H and Hj are re- 
lated by a non-bijective transformation. However, since 
Hj and Hji are in one-to-one correspondence, we can 
write TL = 0^~q Hj = Hk © 14 {dimHk = 00 and 
dimVfe = fc) . Using this notation the multiphoton annihi- 
lation operators are unitarily related to the original anni- 
hilation operator (acting on Hk) via the unitary transfor- 
mation U = J2 n md n ) k ® I m ))(kn + m I [li| performing 
the transformation from H to Hk © In particular, 
UAWW =a©7. 

The Gaussian states are naturally a special case of 
these non-Gaussian states as one would expect. For two 
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FIG. 1: (Color online) The Wigner function of the reduced 
state of a fc-photon two-mode squeezed vacuum state as rep- 
resented in the phase space defined by x\ — a\ + a,pi = 
i[a\ — &i) with k = 3. This state is non-Gaussian. 
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the entanglement. This holds for bjj, but also for mul- 
tiphoton operators A^ k \ hence the criterion derived for 
standard Gaussian states can be directly applied to mul- 
tiphoton Gaussian states as it was demonstrated in the 
present work. A question remains is how to efficiently 
verify whether a given state belongs to a certain sector 
of the Hilbert space TLj for a given k, or not. The answer 
can be given by analyzing the expression for the state un- 
der consideration in the Fock basis. If the populated (i.e. 
non-vanishing) levels are separated by the same energy 
(equivalently, by the same number of photons), then the 
state belongs to a multiphoton sector Ti of the Hilbert 
space H. and its multiphoton Wigner function can be fur- 
ther analyzed. 
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FIG. 2: (Color online) The Wigner function of the reduced 
state of a fc-photon two- mode squeezed vacuum state as repre- 
sented in the phase space defined by X\ — A{ + A^ k \ A = 
i{A[ k)t - A[ k) ) with fe = 3. The state is a thermal state and 
its Gaussian nature is clearly evident. 

modes the operation moving from the basis in which the 
states have a Gaussian Wigner representation to that in 
which they don't is local unitary and as such preserves 



FIG. 3: The maximally entangled mixed states of the multi- 
photon Gaussian states plotted at the same average energy 
(E) — 1 (units are dimensionless) for different values of k. 
The log- negativity, £, is used to measure the entanglement 
while the purity, P, is used to indicate how mixed the state 
is. 
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